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ABSTRACT 


A  new  methodology  has  been  developed  for  applica¬ 
tion  of  Kalman  Filter/Smoothers  to  post-flight  processing  of 
helicopter  flight  test  dynamic  measurements.  This  process¬ 
ing  includes  checking  for  kinematic  compatibility  among  the 
measurements,  identification  of  a  measurement  error  model, 
and  reconstruction  of  both  measured  and  unmeasured  time 
histories.  Emphasis  is  placed  on  identification  of  a  paramet¬ 
ric  measurement  error  model  which  is  valid  for  a  set  of  flight 
test  data.  This  is  facilitated  through  a  new  method  of  con¬ 
catenating  several  maneuver  time  histories.  The  method¬ 
ology  also  includes  a  model  structure  determination  step 
which  ensures  that  a  physically  realistic  parameterization  has 
been  achieved.  Application  of  the  methodology  to  a  set  of 
BO- 105  flight  test  data  is  illustrated.  The  resulting  mini¬ 
mally  parameterized  error  model  is  shown  to  characterize  the 
measurement  errors  of  the  entire  data  set  with  very  little  vari¬ 
ation  in  the  parameter  values.  Reconstructed  time  histories 
are  shown  to  have  increased  bandwidths  and  signal-to-noise 


ratios.  ~ 
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NOMENCLATURE 

°*>  °m 

=  body  axis  translational  accelerations  (m/sec2) 

b 

=  bias  error 

B 

=  tracking  bearing  angle  (rad) 

c.g. 

=  center  of  gravity 

E 

=  tracking  elevation  angle  (rad) 

f,9 

=  vector  functions 

9 

=  local  gravitational  acceleration  (m/sec2) 

P,Q,r 

=  body  axis  angular  velocities  (rad/sec) 

R 

=  tracking  range  (rad) 

t 

=  time 

U,  V,  w 

■  body  axis  translational  velocities  (m/sec) 

u 

=  input  vector 

V 

=  total  velocity  (m/sec) 

Prefer  ted  it  the  46th  Annual  Forum  of  the  American  Helicopter  Society, 
Washington,  DC,  May  1990. 
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=  Earth  axis  displacements  w.r.t.  tracking  station 

(m) 

=  state  vector 
=  observation  vector 
=  angle  of  attack  (rad) 

=  angle  of  sideslip  (rad) 

=  coherence  function 
=  pilot  control  deflection  (%) 

=  random  error 
=  error  parameter 
=  scale  factor 
=  measured  variable 
=  time  delay  (sec) 

=  Euler  angles  (rad) 


INTRODUCTION 


High-quality  flight  test  data  is  required  for  many  appli¬ 
cations  in  rotorcraft  dynamics,  controls,  guidance,  naviga¬ 
tion,  performance,  and  handling  qualities  analyses.  How¬ 
ever,  obtaining  accurate  and  complete  flight  test  data  is  a 
difficult  task.  Measurements  are  always  subject  to  sources 
of  systematic  and  random  error.  In  addition,  some  quanti¬ 
ties  may  be  difficult  to  measure  directly.  This  is  particularly 
true  for  helicopters  because  of  the  their  high  vibration  levels, 
and  is  true  of  VTOL  aircraft  in  general  because  of  their  large 
ranges  of  motion  compared  to  conventional  aircraft. 

Tb  achieve  the  necessary  level  of  quality,  errors  in  the 
flight  test  data  must  be  detected  and  removed.  This  is  pos¬ 
sible  because  the  kinematic  quantities  measured  in  flight  are 
related  to  each  other  by  well-known  differential  equations. 
The  information  in  these  equations  can  be  used  to  determine 
the  kinematic  consistency  (or  compatibility)  of  the  various 
measurements  and  to  reconstruct  others.  This  process  is  of¬ 
ten  called  data  consistency  analysis  and  has  been  used  as 
a  precursor  to  such  activities  as  system  identification  (1,2), 
simulation  validation  (3),  and  accident  investigations  (4). 
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Each  of  these  methods  has  its  own  advantages  and  dis¬ 
advantages.  The  objective  of  all  of  them,  however,  is  to  de¬ 
termine  a  model  of  flight  data  measurement  emus.  The  va¬ 
lidity  of  this  model  is  of  paramount  importance  because  it 
will  affect  the  validity  of  all  subsequent  analyses. 

This  paper  presents  a  comprehensive,  systematic 
methodology  for  conducting  data  consistency  analyses  us¬ 
ing  a  previously  developed  state-estimation  algorithm.  The 
methodology  was  developed  with  the  goal  of  achieving  a  ro¬ 
bust  technique  for  determining  valid  error  models,  hi  the 
process  several  new  techniques  were  developed  which  may 
be  of  interest  to  others  working  in  the  field.  The  application 
of  the  methodology  to  a  set  of  BO-105  flight  test  data  is  il¬ 
lustrated,  and  the  validity  and  usefulness  of  the  results  are 
demonstrated. 

The  flight  test  data  were  produced  by  DLR,  Institute 
for  Flight  Mechanics,  Braunschweig,  Germany,  for  use  by 
AGARD  Flight  Mechanics  Panel  Working  Group  18,  “Ro- 
torcraft  System  Identification,”  December  1987,  and  are 
used  with  permission. 

BACKGROUND 

Kinematic  Equations 

The  set  of  kinematic  equations  for  motions  of  the  air¬ 
craft  center  of  gravity  with  respect  to  a  flat,  non-rotating  earth 
are  presented  in  equations  1  through  IS. 


u  =  rv  —  qw  +  az  -  gsin  6 

(1) 

»  =  pu)-ru+a,+  9  cos  9  sin  ^ 

(2) 

w  =  qu  —  pv  +  a,  +  g  cos  6  cos  <f> 

(3) 

^=p+gsin<£tan0+rcos^tan0 

(4) 

9  -  q  cos  <f>  —  r  sin  <p 

(5) 

=  7  sin  4>sec  9  +  rcos  <£sec  6 

(6) 

X  =  u cos  0 cos  +  v(sin  <f> sin  flcos  -  cos  <j> sin  VO  (7) 


+w(cos^sin0cosV>+  sin  4>  sin  rfi) 

Y  =  ucos*sinV  +  «(sin ^sin 0 sin V>  +  cos^cos (8) 


+ w(cos  4>  sin  6  sin  -  sin  ^  cos  V>) 

Z  =  -usinfl  +  v cos 0 sin  <j>+  wcos  9cos  <t>  (9) 

V  -  W  +  v2  +  w2  (10) 

a  =  arctan(ui/u)  (11) 

P  *  arctan(v/u)  (12) 

R=VX2  +  Y2  +  Z2  (13) 

E  *  arcsin (-Z/R)  (14) 

B  =  arctan(y/JT)  (15) 


For  data  consistency  analyses,  the  equations  are  usually 
considered  in  the  form  of  equations  16  and  17  with  state. 


input,  and  observation  vectors  as  defined  in  equations  18 
through  20. 


i  =  /(x,u) 

(16) 

*  =  p(u) 

(17) 

**■-[«  ,v,w,<j>,9,4>,X,Y,Z) 

(18) 

«*=  [ax,o,,a1,p,9,r] 

(19) 

(20) 

Measurement  Errors 

Measurement  errors  can  be  broadly  classified  as  either 
deterministic,  or  random  errors  as  shown  in  Figure  1.  Bias 
and  scale  factor  errors  generally  result  from  incorrect  in¬ 
strument  calibrations.  Typical  sources  of  time  skew  errors 
are  analog  signal  conditioning  and  multiplexing  of  digital 
data.  Dropouts  are  usually  present  only  in  telemetered  data 
or  tracking  data  Disturbance  errors  are  caused  when  a  sen¬ 
sor  measures  a  motion  not  associated  with  the  rigid-body 
rotation  or  translation  of  the  aircraft  center  of  gravity.  A 
common  example  of  disturbance  errors  is  the  effect  of  atmo¬ 
spheric  turbulence  on  the  measurement  of  the  air-data  vari¬ 
ables.  Quantization  errors  are  caused  by  the  analog-to-digital 
conversion  process  with  finite  digital  word  length.  Other 
sources  of  random  error  include  sensor  and  instrumentation 
noise  and  high  frequency  vibrations. 

Measurement  Error  Models 

A  linear  parametric  model  of  measurement  errors  is 
generally  adopted  for  data  consistency  analysis.  Equations 
such  as 

£«(*)  =  *{  ■  Cm(t  —  **)  +  if  +  *te  (21) 

relate  the  measured  time  histories,  (m,  to  the  estimated  time 
histories,  £„  through  such  parameters  as  bias  errors,  b(,  scale 
factors,  time  delays,  t{,  and  random  errors,  r%.  Some 
methods  do  not  model  the  random  error  explicitly,  and  those 
that  do  generally  assume  it  to  be  white  gaussian.  Time  delay 
errors  are  seldom  modeled,  but  some  researchers  have  real¬ 
ized  significant  improvements  in  results  by  including  them 
in  the  error  model  (9,10). 

Data  Consistency  Analysis 

Parameter  values  are  determined  by  fitting  the  measured 
time  histories  to  the  estimated  time  histories  subject  to  the 
constraints  of  equations  16  and  17.  Differential  equations 
for  the  deterministic  error  parameters  are  appended  to  to  the 
state  equations  and  provide  additional  constraints.  These  are 
usually  trivial  since  the  error  parameters  are  usually  assumed 
to  be  time-invariant. 

Algorithms-  Several  algorithms  have  been  used  to 
identify  parameter  values  and  calculate  estimated  time  histo¬ 
ries.  They  differ  mainly  in  the  level  of  sophistication  of  their 
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Fig.  1  Types  of  measurement  errors. 


constraints  of  equations  16  and  17.  Differential  equations 
for  the  deterministic  error  parameters  are  appended  to  to  die 
state  equations  and  provide  additional  constraints.  These  are 
usually  trivial  since  the  error  parameters  are  usually  assumed 
to  be  time-invariant. 

Algorithms-  Several  algorithms  have  been  used  to 
identify  parameter  values  and  calculate  estimated  time  histo¬ 
ries.  They  differ  mainly  in  the  level  of  sophistication  of  their 
random  error  models  and  the  type  of  cost  function  they  em¬ 
ploy.  Some  of  the  more  common  algorithms  are  discussed 
below. 

Least  Squares-  Random  errors  are  not  modeled  ex¬ 
plicitly  in  this  method.  The  state  equations  are  integrated  us¬ 
ing  measured  state  and  input  variables  to  calculate  estimates 
of  the  state  and  output  variables.  This  process  continues,  us¬ 
ing  estimated  state  values  in  further  iterations,  until  param¬ 
eter  values  are  determined  which  minimize  the  least  squares 
fit  error  between  the  measured  and  estimated  outputs. 

Extended  Kalman  Filter-  Random  errors  in  both  the 
input  variables  (process  noise)  and  output  variables  (mea¬ 
surement  noise)  can  be  modeled  explicitly  in  this  method. 
The  presence  of  process  noise  means  that  the  state  variables 
must  be  treated  as  stochastic  quantities.  The  Kalman  filter 
steps  through  the  data  calculating  estimates  of  the  state  and 


measurement  variables  based  on  models  of  the  system  dy¬ 
namics  and  statistical  properties  of  the  process  and  measure¬ 
ment  noise.  Parameter  values  are  determined  through  mini¬ 
mization  of  a  cost  function  based  on  state  and  measurement 
variable  fit  errors  weighted  by  the  process  and  measurement 
error  covariance  matrices  respectively. 

Maximum  Likelihood-  The  maximum  likelihood 
identification  technique  is  usually  applied  to  the  data  consis¬ 
tency  analysis  problem  as  an  output  error  (process  noise  ne¬ 
glected)  formulation  (8,9,10).  This  allows  the  state  equations 
to  be  integrated  directly.  Random  measurement  errors  are 
modeled  explicitly  by  a  covariance  matrix.  Parameter  val¬ 
ues  are  determined  through  minimization  of  the  likelihood 
cost  function  which  includes  a  weighted  least  squares  mea¬ 
surement  fit  error  term  and  an  additional  logarithmic  term. 

Model  Structure  Determination 

The  model  structure  refers  to  the  set  of  non-zero  pa¬ 
rameters  which  are  used  to  parameterize  the  measurement 
error.  Determination  of  an  appropriate  model  structure  is 
an  important  part  of  any  parameter  identification  problem. 
Improvements  in  the  fit  can  almost  always  be  achieved  by 
adding  non-zero  parameters  to  the  model.  However,  only 
parameters  which  truly  represent  a  physical  characteristic  of 
the  system  should  be  included.  An  incorrect  model  structure 
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will  cause  the  parameter  estimates  to  be  biased  away  from 
their  true  values. 

TWo  factors  which  should  be  considered  during 
model  structure  determination  are  parameter  insensitivity 
and  parameter  correlation.  If  the  solution  (cost  function)  is 
relatively  insensitive  to  variations  in  a  particular  parameter, 
then  this  parameter  is  said  to  have  a  high  insensitivity.  The 
parameter  is  then  less  likely  to  be  an  important  part  of  the 
model  structure,  and  removing  it  should  be  considered 

If  variations  in  two  or  more  parameters  cause  sim¬ 
ilar  variations  in  the  solution  (cost  function),  then  they  are 
said  to  be  correlated  The  more  two  parameters  are  corre¬ 
lated,  the  less  possible  it  is  to  determine  unique  values  for 
each  of  them.  This  in  turn  makes  it  difficult  to  determine 
whether  each  of  them  is  an  important  part  of  the  model.  This 
can  be  remedied  by  fixing  one  of  them  at  an  a  priori  value 
if  enough  information  is  available  to  determine  a  reasonable 
one. 

METHODOLOGY 

A  comprehensive  methodology  has  been  developed 
for  application  of  Kalman  Filter/Smoothers  to  the  identifica¬ 
tion  of  models  for  helicopter  flight  test  data  measurement  er¬ 
rors.  A  flow  chart  of  the  information  flow  through  this  sys¬ 
tem  is  shown  in  Figure  2.  The  filter/smoother  algorithm  is 
applied  in  both  the  initial  solution  and  final  solution  blocks. 
Some  of  the  major  features  of  the  methodology  are: 

1.  All  of  the  flight  test  data  for  a  given  flight  test  is 
passed  through  a  preprocessor  whose  function  is  to  increase 
the  signal-to-noise  ratio  of  the  data  set. 

2.  A  small  subset  of  the  preprocessed  data  is  processed 
in  an  initial  solution  to  obtain  an  initial  error  model  structure 
for  use  as  a  start-up  in  the  final  solution.  Knowledge  of  in¬ 
strumentation  characteristics  is  used  at  this  step  to  determine 
appropriate  values  for  the  measurement  error  covariance  ma¬ 
trix  of  the  Kalman  Filter. 

3.  All  of  the  flight  test  data  is  then  processed  in  the  final 
solution  step  to  obtain  an  error  model  valid  for  the  entire  set 
of  data. 

The  final  error-model  parameters  can  be  applied  directly 
as  corrections  to  the  raw  flight  test  data,  or  the  filter/smoother 
can  be  used  with  the  error  parameters  fixed  to  these  values 
to  reconstruct  estimates  of  error-free  time  histories. 

Each  of  the  major  sections  of  this  process  are  explained 
in  greater  detail  in  the  sections  which  follow. 

SMACK 

For  the  present  study,  the  filter/smoother  program 
SMACK  (SMoothing  for  AirCraft  Kinematics),  developed 
by  Ralph  Bach  at  Ames  Research  Center,  was  employed  as 
the  parameter  estimation/state  reconstruction  algorithm  in 
the  initial  and  final  solution  blocks.  A  block  diagram  of  the 


Fig.  2  Flow  chart  of  developed  methodology. 


state  and  measurement  models  used  by  SMACK  is  shown 
in  Figure  3.  These  models  differ  from  the  more  commonly 
employed  models  of  equations  18  through  20  because  they 
have  been  formulated  to  keep  the  state  model  linear.  All  mea¬ 
sured  variables  are  also  treated  as  observations  and  never  as 
system  inputs.  Estimation  of  constant  bias  and  scale  factor 
errors  in  each  of  the  measurements  as  well  as  time-varying 
winds  is  possible.  Other  algorithms  could  be  substituted  for 
SMACK  without  reducing  many  of  the  benefits  of  the  over¬ 
all  methodology.  However,  SMACK  has  several  properties 
and  features  which  n-A iu.  >t  a  favorable  choice.  These  are 
discussed  below. 

Stochastic  Optimal  Estimator-  The  SMACK  algo¬ 
rithm  consists  of  a  backward  “information”  filter  and  forward 
smoother.  This  formulation  offers  all  of  the  advantages  of 
the  extended  Kalman  filter  plus  certain  computational  advan¬ 
tages  which  result  in  a  more  accurate  and  less  computation¬ 
ally  expensive  solution  (6).  The  use  of  a  linear  state  model 
also  improves  the  performance  of  the  SMACK  algorithm. 
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Fig.  3  State  and  measurement  models  used  by  SMACK. 
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Other  Computational  Features-  SMACK  calculates 
the  means  and  standard  errors  of  the  final  measurement  resid¬ 
uals,  which  supplies  a  useful  feedback  to  the  user  about  the 
accuracy  of  initial  estimates  of  the  measurement  error  stan¬ 
dard  deviations.  A  provision  is  also  supplied  for  the  elimi¬ 
nation  of  selected  data  points  from  the  measurement  update. 
This  is  very  useful  for  elimination  of  disturbance-type  errors 
from  the  measurements. 

Convenience  Features-  A  predefined  measurement 
model  is  available,  but  new  measurements  are  easily  added. 
An  internal  simulated  maneuver  is  available  for  debugging 
new  input  decks.  Printed  output  provides  many  useful  met¬ 
rics  of  solution  quality;  and  graphical  output  of  measure¬ 
ment,  estimate,  and  residual  time  histories  is  supplied. 

Preprocessing 

Studies  have  been  conducted  to  investigate  the  variabil¬ 
ity  of  parameter  values  identified  from  different  flight  test 
maneuvers  (1,7).  This  variability  is  often  caused  by  poor  ex¬ 
citation  of  some  states  for  certain  types  of  maneuvers.  Klein 
(8)  postulates  that  it  may  be  possible  to  increase  the  accu¬ 
racy  of  parameter  estimates  by  designing  an  optimal  maneu¬ 
ver.  Indeed,  this  approach  has  been  taken  in  a  flight  test  pro¬ 
gram  of  the  AV-8B  Harrier  aircraft  conducted  at  the  Ames 
Research  Center  (2).  However,  it  is  often  the  case  that  flight 
tests  are  conducted  without  this  consideration  in  mind.  The 
approach  taken  here  is  to  create  “optimal"  maneuvers  from 
sub-optimal  flight  test  records  in  a  preprocessing  step. 

Flight  test  time  histories  from  several  different  types  of 
maneuvers,  usually  one  for  each  control  axis,  are  concate¬ 
nated  to  construct  a  composite  maneuver  with  sufficient  ex¬ 
citation  of  all  kinematic  quantities.  Discontinuities  at  the 
junctions  between  events  are  avoided  by  inserting  a  buffer  re¬ 
gion  of  time  between  the  records.  These  sections  are  flagged 
for  elimination  from  the  filter  measurement  update  allow¬ 
ing  SMACK  to  construct  trajectories  joining  the  two  events. 
These  “pseudo-dropouts”  are  of  sufficient  length  to  allow  the 
estimated  aircraft  states  to  transition  smoothly.  When  track¬ 
ing  data  are  available,  coordinate  transformations  are  made 
to  the  data  to  ensure  a  continuous  flight  path  throughout  the 
composite  maneuver. 

All  of  the  flight  test  data  are  also  carefully  low-pass  fil¬ 
tered  to  remove  noise.  This  approach  is  taken  because  dif¬ 
ferent  types  of  sensors  have  different  bandwidths.  For  in¬ 
stance,  accelerometers  will  be  much  more  sensitive  to  high- 
frequency  vibrations  than  gyros  will  be.  Thus,  low-pass 
filtering  removes  another  possible  source  of  inconsistency 
among  measurements. 

Also,  a  visual  inspection  of  the  data  is  made  and  dropout 
and  disturbance  errors  are  flagged  for  elimination  from  the 
filter  measurement  update.  This  has  the  effect  of  increasing 
the  signal-to-noise  ratio  of  each  affected  measurement  and 
removing  a  source  of  error  which  could  otherwise  bias  esti¬ 
mates  of  the  error  model  parameters. 


Use  of  this  preprocessing  procedure  has  resulted  in  less 
variation  in  identified  parameter  values  and  improved  time 
history  estimates. 

Initial  Solution 

The  objective  of  the  initial  solution  step  is  to  determine 
an  initial  error  model  structure  and  refine  estimates  of  the 
measurement  error  covariance  values.  This  analysis  is  per¬ 
formed  by  iteratively  processing  one  representative  compos¬ 
ite  maneuver  through  SMACK.  The  results  are  used  as  a 
start-up  for  processing  of  all  of  the  cases  in  the  final  solu¬ 
tion. 

The  initial  values  of  the  measurement  error  covariances 
are  determined  from  the  digital  resolution  of  the  instrumenta¬ 
tion  system.  This  information  is  usually  available  from  cal¬ 
ibration  data  for  each  measured  channel  as  the  number  of 
engineering  units  per  “count." 

For  example,  a  typical  12-bit  system  would  provide  a 
maximum  of  2 12  =  4096  counts  per  channel.  If  an  ac¬ 
celerometer  was  calibrated  to  read  -1 .00g  at  1600  counts  and 
+1.50g  at  2800  counts,  the  digital  resolution  of  that  channel 
would  be  2.50g(1200  *  0.00208g. 

These  values  are  input  to  SMACK  as  the  standard  de¬ 
viations  of  the  random  error  component  of  the  measure¬ 
ment  errors.  This  is  compatible  with  the  requirements  of  the 
filter/smoother  since  quantization  error  can  be  modeled  as 
white  noise  (12,  pages  190-195) .  The  assumption  that  all  of 
the  random  error  may  be  attributed  to  quantization  effects  is 
sufficiently  accurate  for  starting  up  the  initial  solution  pro¬ 
cess.  In  fact,  it  is  likely  that  low-pass  filtering  has  removed 
most  other  sources. 

The  representative  case  is  then  processed  iteratively 
through  SMACK.  After  each  iteration,  the  error  model  struc¬ 
ture  and  error  covariance  values  may  be  modified  based  on 
evaluation  of  the  three  following  criteria. 

Measurement  Residuals-  In  his  forthcoming  NASA 
Reference  Publication  “State  Estimation  Applications  in  Air¬ 
craft  Flight-Data  Analysis  (A  User’s  Guide  for  SMACK),” 
Bach  states  that  one  of  the  criteria  for  a  “good”  SMACK  so¬ 
lution  is  that  all  of  the  residuals  should  be  white,  zero  mean, 
and  have  standard  deviations  close  to  the  measurement  error 
standard  deviation  values  provided  by  the  user.  This  concept 
is  similar  to  the  Innovations  Property  of  linear-optimal  filters 
(13,  page  400) . 

At  each  iteration  of  the  initial  solution  the  residuals 
are  examined  for  whiteness,  offset,  and  magnitude.  Non¬ 
white  or  non-zero  mean  residuals  may  be  indicative  of  error- 
parameters  which  are  missing  from  the  model.  A  missing 
scale  factor  will  cause  the  corresponding  residual  signal  to 
have  a  shape  similar  to  the  measured  signal.  A  missing  bias 
parameter  will  cause  the  residual  to  have  a  non-zero  mean 
value. 
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The  standard  errors  of  the  residuals  are  compared  to  the 
user-supplied  estimates  of  the  measurement  error  standard 
deviations.  If  a  residual  is  sufficiently  white,  then  its  stan¬ 
dard  error  will  be  a  good  measure  of  the  true  measurement 
error  standard  deviation.  If  not,  an  estimate  of  the  stan¬ 
dard  error  of  the  random  portion  of  the  residual  signal  may 
have  to  be  made.  The  user  may  then  update  his  estimates 
accordingly. 

Parameter  Insensitivity-  The  insensitivity  of  the 
model  to  variations  in  its  parameters  is  determined  through 
examination  of  the  parameter  Cramer-Rao  lower  bounds 
generated  by  SMACK.  It  is  well  known  that  the  Cramer-Rao 
bounds  overestimate  the  accuracy  of  the  identified  parame¬ 
ter,  but  they  are  useful  for  comparisons  of  relative  accuracy 
between  parameters,  and  variations  from  case  to  case  (10). 
When  the  ratio  of  the  Cramer-Rao  bound  to  the  parameter 
value  is  relatively  large  (usually  by  a  factor  of  ten)  for  a  par¬ 
ticular  parameter,  it  is  eliminated  from  the  model  structure. 
It  has  been  the  author’s  experience  that  there  is  little  effect  on 
the  whiteness  of  the  residuals  from  elimination  of  insensitive 
parameters. 

Of  course,  physical  understanding  of  a  particular  prob¬ 
lem  should  always  be  used  to  pare  down  a  model  structure 
from  the  outset.  For  example,  if  one  is  conducting  an  anal¬ 
ysis  of  the  three  rotational  degrees-of-freedom  only,  the  so¬ 
lution  will  be  very  insensitive  to  bias  errors  in  ^  and  0.  For 
a  six-degree -of-freedom  analysis  without  tracking  data  the 
solution  will  be  completely  insensitive  to  a  bias  error  in  ip. 
This  is  readily  verified  through  inspection  of  the  kinematic 
equations. 

Parameter  Correlation-  No  numerical  indicator  of  pa¬ 
rameter  correlation  is  presently  calculated  in  SMACK.  Cor¬ 
relation  of  parameters  is  currently  avoided  by  consideration 
of  the  the  structure  of  the  kinematic  equations  and  use  of  non- 
kinematic  information.  For  example,  bias  errors  in  a*  and  0 
will  be  highly  correlated.  The  kinematic  equations  reveal 
the  correlation,  but  more  information  is  needed  to  determine 
which  measurement  is  really  biased.  If  the  flight  test  data 
includes  segments  of  trimmed  flight,  then  pilot  comments  or 
simulation  results  may  be  used  to  determine  the  validity  of 
the  measured  trim  attitude. 

Convergence-  The  above  criteria  are  examined  at  each 
iteration  until  an  appropriate  model  structure  and  measure¬ 
ment  error  covariance  values  are  obtained.  Convergence 
problems  have  not  been  encountered. 

Final  Solution 

The  direction  in  which  the  analysis  proceeds  from  this 
point  depends  on  the  characteristics  of  the  database.  If  the 
flight  tests  were  conducted  over  a  long  period  of  time,  then  it 
is  likely  that  the  instrumentation  characteristics  will  vary  sig¬ 
nificantly  throughout  the  database.  In  this  case  we  can  expect 
the  error  model  structure  and  parameter  values  to  vary,  and 
the  procedure  of  the  initial  solution  will  have  to  be  applied 


enough  times  to  capture  this  variation.  However,  if  the  flight 
tests  are  conducted  within  a  short  enough  period  of  time,  the 
assumption  of  a  stationary  error  model  can  be  made.  In  this 
case  the  procedure  of  the  final  solution  is  employed  to  de- 
termine  a  model  structure  and  parameter  values  valid  for  the 
entire  data  base. 

In  the  final  solution  step,  SMACK  is  applied  iteratively 
to  all  of  the  composite  maneuver  cases.  At  each  iteration,  the 
error-model  structure  is  held  constant  across  all  the  maneu¬ 
vers.  The  model  structure  from  the  initial  solution  is  used  as 
a  start-up  in  this  process. 

After  each  iteration,  the  sample  mean  values  and  stan¬ 
dard  deviations  of  the  parameters  are  calculated.  The  pa¬ 
rameter  with  the  largest  variation  above  a  certain  threshold 
is  dropped  from  the  model.  This  process  continues  until  all 
of  the  parameters  have  sufficiently  small  variations. 

It  has  been  the  author’s  experience  that  an  appropriate 
threshold  value  of  parameter  variation  usually  appears  dur¬ 
ing  the  analysis.  The  parameters  are  usually  grouped  into 
sets  with  low-variance  and  high-variance.  Sequential  elim¬ 
ination  of  the  high-variance  parameters  tends  to  reduce  the 
variations  in  the  low-variance  parameters,  accentuating  this 
division.  Determination  of  an  appropriate  threshold  value  is 
therefore  a  simple  matter. 

When  a  final  solution  is  determined,  the  model  struc¬ 
ture  and  parameter  values  are  enforced  in  SMACK  to  recon¬ 
struct  time  histories  of  measured  and  unmeasured  variables 
for  each  of  the  unlinked  flight  test  events. 

RESULTS 

The  proposed  methodology  has  been  applied  to  a  set 
of  flight  test  data  produced  by  the  DLR,  Institute  for  Flight 
Mechanics  in  Braunschweig,  Germany,  for  use  in  time-  and 
frequency-domain  parameter  identification  (1).  These  tests 
were  flown  in  the  DLR  ATTHeS  BO-105  helicopter  (see 
Figure  4)  The  data  include  doublet,  multi-step  “3211", 
and  frequency-sweep  maneuvers  in  each  of  the  four  control 
axes.  The  doublets  were  repeated  twice  in  each  direction;  the 
3211’s  were  repeated  three  times  in  each  direction;  and  the 
frequency-sweeps  were  repeated  three  times  for  a  total  of  52 
flight  test  events.  The  nominal  flight  condition  was  horizon¬ 
tal  flight  at  80  knots  and  a  density  altitude  of  3000  feet,  and 
the  entire  flight  test  was  conducted  in  less  than  two  weeks. 
The  measurements  used  in  consistency  checks  and  state  re¬ 
construction  included:  the  Euler  angles  <p  and  6\  body  angu¬ 
lar  rates  p,  q,  and  r,  c.g.  specific  forces  a*.  a,,  and  at;  and 
c.g.  translational  velocities  u,  v,  and  w.  The  translational 
velocities  had  previously  been  derived  from  measurements 
taken  with  a  HADS  air  data  system  which  uses  a  swiveling 
pitot  static  mechanism  located  below  the  main  rotor. 

The  52  flight  test  events  were  preprocessed  into  1 3  com¬ 
posite  maneuvers,  each  consisting  of  a  longitudinal  cyclic 
maneuver,  a  lateral  cyclic  maneuver,  a  collective  maneuver. 


7 


Fig.  4  DLR  ATTHes  BO105  helicopter. 


and  a  pedals  maneuver.  Disturbance  errors  in  the  velocity 
signals  were  discovered  in  many  of  maneuvers.  It  is  believed 
that  these  were  caused  by  interaction  of  the  HADS  air  data 
sensor  with  the  rotor  wake.  They  were  flagged  for  elimina¬ 
tion  from  the  filter  measurement  update. 

Parameter  values  for  the  final  model  structure  are  tab¬ 
ulated  for  the  13  composite  maneuvers  in  Table  1.  The 
SMACK  calculated  Cramer-Rao  lower-bounds  are  also  tab¬ 
ulated.  The  standard  deviations  of  all  of  the  parameters  are 
less  than  11%  of  their  respective  mean  values.  Such  small 
variations  in  the  parameter  estimates  validate  the  assump¬ 
tion  of  a  stationary  error  model  and  provide  a  high  degree  of 
confidence  in  the  model  structure  and  parameter  mean  values 
themselves.  It  is  also  notable  that  the  Cramer-Rao  bounds 
underestimate  the  sample  variations  by  as  much  as  two  or¬ 
ders  of  magnitude. 

The  usefulness  of  the  composite  maneuver  approach  can 
be  seen  by  comparing  the  results  in  Ttble  1  with  those  in  Ta¬ 
ble  2.  The  parameter  values  in  Table  2  were  obtained  using 
the  same  model  structure,  but  the  individual  flight  test  ma¬ 
neuvers  were  used  in  the  identification  instead  of  the  com¬ 
posite  maneuvers.  It  is  evident  that  the  standard  deviations 


of  the  parameters  increase  and  their  mean  values  change  sig¬ 
nificantly  when  this  approach  is  taken. 

The  effect  of  excitation  of  a  particular  variable  on  the 
identification  of  error  parameters  associated  with  it  can  be 
seen  in  Table  3.  Here  the  identified  values  for  the  scale  factor 
on  lateral  velocity  are  tabulated  by  the  control  axis  of  the 
maneuver.  It  is  clear  that  variations  in  the  identified  value 
of  X,  are  smaller  for  the  pedal  and  lateral  cyclic  maneuvers 
which  contain  much  larger  lateral  velocity  responses  than  the 
other  maneuvers. 

Comparison  of  filtered  flight  test  data  and  reconstructed 
time  histories  for  a  typical  composite  maneuver  is  illus¬ 
trated  in  Figure  5.  7116  reconstructed  data  have  been  de- 
corrected  by  the  identified  error  parameters  for  plotting  pur¬ 
poses.  Also,  the  “pseudo-dropouts”  between  events  have 
been  removed  to  improve  the  plot  scales.  Only  lateral  vari¬ 
ables  are  plotted  to  save  space,  but  the  results  are  represen¬ 
tative  of  the  the  fits  in  other  axes. 

The  consistency  is  generally  very  good.  The  lack  of  re¬ 
maining  deterministic  errors  indicates  that  the  error  model 
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Table  1  Error  mode!  parameters  identified  for  composite  maneuvers 


■KTC/DBHi 

I - 

mm  ii  imiiiB 

Hmucisaurcctiw 

mvriMtvAwsm 

L 

■WIL&iPattUJJM 

■wjr£'<aaiiM:m 

■iXrMt'.xaiXMIlH 

wA:i*iFa«:uw 

3 

■UMraaixuiiM 

4 

wxiswxaiMwmi 

■IlMHCaiXMIlH 

«x:i2ii:aixw)MH 

"  5 

mwxuzuMmm 

■*jt?:.r»ea»xi:i:»x«i 

■ir.twetirm/* 

6 

mWMfiVAU&M 1 

wxittEsaixiicuM 

WX'H’U’IXOJM 

Mtx.Y£i:sa»x».uiMi 

mxmfVKWim 

7 

mmz&wimm 

■«t»,*isaw;i;i« 

wjoea*XM:i»M 

mwsnsanKmm 

my&mxvALLLim 

a 

MamiEwmm 

»aw«txi:t:<)r*B 

smmiiiamrxMm 

Mtx:iy>;aix*:w 

y 

Muwiv&iatMtm 

mam&viLLitm 

muisM&aimm 

■IHWiaillllM 

10 

■qWMKWW 

■wwssmw 

■umishuumm 

■ttaiatgaaaiM 

■  li  1 

mumsasamim 

17 

■uxiKwsawiw 

MLiSi/Aiix 

■imubbimm 

13 

■mwnww:w 

mmrr.uszmjitm 

mam&Hasaum 

0.9556*0.0272 

0.8750*0.0378  1 

adequately  characterizes  those  present  in  the  flight  data  Dis¬ 
turbance  errors  are  evident  in  the  lateral  velocity  measure¬ 
ment.  Vertical  dotted  lines  indicate  the  regions  of  the  sig¬ 
nal  which  have  been  eliminated  from  the  measurement  up¬ 
date.  The  missing  velocity  segments  have  been  satisfactorily 
reconstructed. 

Improvements  in  the  velocity  signals  are  also  clearly 
seen  in  the  frequency  domain.  The  frequency  responses  for 
lateral  velocity  due  to  pedal  inputs  identified  using  measured 
and  reconstructed  lateral  velocity  signals  are  compared  in 
Figure  6.  Bode  magnitude  and  phase  and  the  coherence  func¬ 
tion  72  are  plotted.  The  coherence  function  indicates  the 
fraction  of  the  output  which  can  be  accounted  for  by  linear 
relation  with  the  input  It  may  be  reduced  from  a  maximum 
value  of  1  by  nonlinearities,  input  and  output  noise,  and  lack 
of  excitation  and/or  response.  The  most  striking  feature  of 
Figure  6  is  that  the  coherence  function  for  the  response  iden¬ 
tified  from  the  reconstructed  data  is  larger  over  a  much  wider 
frequency  range.  This  is  the  result  of  two  factors.  The  first  is 
that  the  signal-to-noise  ratio  of  the  lateral  velocity  signal  has 
been  increased  through  the  elimination  of  random  and  distur¬ 
bance  errors.  The  second  is  that  the  bandwidth  of  the  lateral 
velocity  signal  has  been  increased  by  the  complementary  fil¬ 
ter  action  of  the  Kalman  Filter.  This  results  from  the  use  of 
many  measurements  with  different  response  characteristic  in 
the  calculation  of  the  reconstructed  time  histories. 

It  is  also  interesting  to  note  that,  for  the  range  of  fre¬ 
quencies  where  both  frequency  responses  have  high  coher¬ 
ence  values,  their  response  characteristics  are  very  similar. 
The  only  major  difference  is  a  magnitude  shift  of  approx¬ 
imately  3dB  which  agrees  with  the  A,  =  0 .71  identified  in 
the  time  domain.  This  also  indicates  that  the  state-estimation 
process  has  not  corrupted  the  dynamic  relationship  between 
the  control  and  response  variables. 

The  results  presented  here  have  been  used  in  time- 
and  frequency-domain  identification  and  verification  of 
BO-105  dynamic  models.  Improvements  in  the  identifica¬ 
tions  wete  obtained  for  use  of  the  identified  parameters  in 
the  time-domain  method  and  use  of  reconstructed  data  in  the 
frequency-domain  method  (1). 


Measured  data 
Reconstructed  data 


Fig.  6  Comparison  of  frequency  responses  identified  from 
flight  test  and  reconstructed  data. 

CONCLUDING  REMARKS 

A  comprehensive  and  systematic  methodology  for  iden¬ 
tifying  models  of  helicopter  flight  data  measurement  errors 
has  been  developed.  This  methodology  makes  use  of  the  fa¬ 
vorable  properties  and  features  of  the  optimal  filter/smoother 
program  SMACK.  A  preprocessing  step  is  used  to  construct 
composite  maneuvers  with  increased  kinematic  information 
content  and  to  increase  the  signal-to-noise  ratios  of  individ¬ 
ual  measurements.  Appropriate  measurement  error  covari¬ 
ance  values  are  determined  in  an  initial  solution  step,  and  a 
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Table  2  Error  model  parameters  identified  for  flight  test 
maneuvers 


Case 

ip  (d/s) 

i,  (d/s) 

A. 

A, 

|— 

1 

0.0933 

0.0930 

0.9472 

0.6684 

0.8181 

2 

0.0835 

0.1016 

1.0010 

0.9543 

0.8986 

*7 

0.0895 

0.0814 

0.9465 

0.6760 

0.7989 

4 

0.0825 

0.0697 

0.9880 

1.0080 

0.9206 

5 

0.0353 

0.0691 

1.0390 

0.6840 

0.8834 

6 

0.0414 

0.0959 

1.1030 

0.8070 

1.0140 

7 

0.0920 

0.0783 

0.9392 

0.8276 

0.8185 

8 

0.1054 

0.1025 

1.0100 

0.7776 

0.9546 

9 

0.0803 

0.0817 

0.9969 

1.0610 

0.8723 

10 

0.0719 

0.0771 

0.8741 

0.3801 

0.6621 

11 

0.0607 

0.0734 

0.9432 

0.7079 

0.8363 

12 

0.1087 

0.0873 

0.9648 

0.5633 

0.7920 

13 

0.0890 

0.0853 

0.9511 

0.6546 

0.8519 

14 

0.1278 

0.0753 

1.1560 

0.7695 

0.9907 

15 

0.0625 

0.1246 

1.0140 

0.7400 

0.8646 

16 

0.0562 

0.0832 

1.0610 

0.7548 

0.8791 

17 

0.1717 

0.2329 

1.0230 

1.6270 

0.8639 

18 

0.0955 

0.0798 

0.9454 

0.6182 

0.8548 

19 

0.1454 

0.0792 

0.9460 

0.9846 

0.8053 

20 

0.1115 

0.1292 

0.8733 

0.3846 

0.7930 

21 

0.0692 

0.0875 

0.8815 

0.9083 

0.7627 

22 

0.0860 

0.1099 

0.8918 

0.5254 

0.8839 

23 

0.0839 

0.0825 

0.9895 

0.5860 

0.8744 

24 

0.0922 

0.1016 

0.9693 

0.6867 

0.9216 

25 

0.0852 

0.1214 

0.9390 

0.4398 

0.8797 

26 

0.1142 

0.1572 

0.9745 

0.7874 

0.8078 

27 

0.0753 

0.1118 

0.9119 

0.9411 

0.9201 

28 

0.0791 

0.1448 

0.9119 

0.8493 

0.9372 

29 

0.0896 

0.1132 

0.9660 

0.8372 

0.9252 

30 

0.0789 

0.0779 

0.9132 

0.9084 

0.8790 

31 

0.0752 

0.0828 

0.9005 

0.8295 

32 

0.0821 

0.0757 

0.9601 

0.7037 

0.9101 

33 

0.0834 

0.1083 

0.9151 

0.5536 

0.8495 

34 

0.0838 

0.0816 

0.9564 

0.6686 

0.7682 

35 

0.0897 

0.0760 

0.9436 

0.6310 

0.8044 

36 

0.1026 

0.1092 

0.9276 

0.6411 

0.8165 

37 

0.0676 

0.1209 

1.0010 

0.6540 

0.8345 

38 

0.1123 

0.0939 

0.9727 

0.8419 

1.0580 

39 

0.0872 

0.1079 

0.9306 

0.5956 

0.7921 

40 

0.0894 

0.0908 

0.9616 

0.6751 

0.9224 

41 

0.0569 

0.0930 

0.8370 

0.5546 

0.7126 

42 

0.0756 

0.0928 

0.9479 

0.7819 

0.8639 

43 

0.0750 

0.0933 

0.8803 

0.6354 

0.7844 

44 

0.0745 

0.0933 

0.7987 

0.4792 

0.7732 

45 

0.0692 

0.0925 

0.8289 

0.5846 

0.8208 

46 

0.0742 

0.0918 

0.9143 

0.6473 

0.8498 

47 

0.0961 

0.0909 

0.9895 

0.7911 

1.1930 

48 

0.1069 

0.0950 

1.0360 

1.0310 

1.1300 

49 

0.0843 

0.0918 

0.8794 

0.6546 

0.7299 

50 

0.0799 

0.0904 

0.8518 

0.5467 

0.6609 

51 

0.0728 

0.0913 

0.9640 

0.7113 

0.7628 

52 

0.1044 

0.0940 

1.0880 

0.7037 

0.8520 

6 

0.0849 

0.0955 

0.9488 

0.6764 

0.8432 

as 

0.0230 

0.0265 

0.0687 

0.2095 

0.1031 

Table  3  Variation  in  A„  identification  with  level  of  lateral 
velocity  excitation  for  different  control  axis  inputs 
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measurement  error  model  valid  for  a  set  of  flight  test  data  is 
identified  in  a  final  solution  step. 

This  methodology  had  been  applied  to  the  identification 
of  a  physically  realistic  measurement  error  model  for  a  set  of 
BO- 105  flight  test  data.  The  identified  model  is  described  by 
a  small  number  of  parameters  which  each  vary  less  than  1 1  % 
over  the  entire  data  base.  The  use  of  composite  maneuvers 
in  the  identification  was  shown  to  decrease  variations  in  the 
parameter  values. 

Reconstructed  time  histories  were  shown  to  be  free  of 
many  deterministic  and  random  errors.  This  included  the 
elimination  of  rotor  wake  interference  effects  from  the  ve¬ 
locity  measurements.  Frequency  response  analysis  demon¬ 
strated  that  these  signals  have  increased  bandwidths  and 
signal-to-noise  ratios. 

Use  of  the  results  of  the  BO-105  data  consistency  analy¬ 
sis  in  identification  of  BO- 105  dynamic  models  has  produced 
improvements  in  the  identification  results. 
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